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A popular choice for modelling turbulence subgrid tensor consists of intro- 
ducing the concept of "eddy viscosity" vturb- 

Tij — -SijTkk = —2vturbSij (1) 

where Sij is the shear stress filtered on the scale A. The quantity vturb characterizes 
the turbulence. Its order of magnitude is given by the product of a length times a 
velocity characteristic of the turbulent motion: vturb '-^ I • Vi. Choosing the length 
equal to the computational grid, / — A, and Va = jS"] A, we obtain the Smagorinsky 
model [1]: 

uturb = CA^ (2) 
The quantity C can be calculated dynamically from the Germano identity: 

which furnishes: ^ 

{SijSij) 

The behaviour of the parameterization ^ is very similar to the dynamic Smagorin- 
sky model formulated in 3 . Hence the tensor given by equation (|4]) is: 



\5^,Tkk = ^-kA^S,, (5) 

{SijSij) 



In order to eliminate the constrain that Tij is parallel to Sij , we propose the following 
parameterization: 

r,, =-2A4/3(^C,j+c^y) (6) 

where is a quantity representative of the strain-rate (for instance one invariant 
of the tensor Sij). The tensor components are calculated from Germano's 
identity directly, without introducing the concept of eddy viscosity, wheras the 
scalar quantity c is the one that better reproduce the dissipation energy. 
Combining ^ and ([3]) we obtain: 

- 2(2A)4/3 cg'A) ^ ^ ^ 2 A4/3 (SC,^ + c%,j ) = (7) 
Assuming that c|J'^^ = c[f^ and Sdj ~ Sdj, we have: 

- 2 A4/3 (^4'^^ + c%,j ) = aUj (8) 
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- ^a^ + 2A*/3c5i (9) 

S S 

In order to calculate c we multiply eq.® by Sij and assume that the following 
relationship is valid for the average quantities: 

(3CA) ^ c{%%,) (10) 
Using the above relationship we have: 

2A4/3c=-iM^ (11) 



^ {Sij Sij) 



Hence the tensor Tij becomes: 



5" a {Li j Sij) f-^ S^ 

i^i^ij + 2 — — — 
S {Sij Sij ) 



iij — a-^Lij + ——^ — ^ — I Sij — ^Sij ) (12) 



and assuming S — \/ Sij Sij we obtain the final form: 



SijSij ^ a{LijSii) (— ySijSij^ 



^ij + o ^ ^ ^ij ; ^ ^ ^ij {^^) 



Sij Sij {Sij Sij) Y \J Sij Sij 



We can see that the tensor contains the first term that is similar to the tensor 
predicted by similarity models [S]: 

'^ij — C sira-^ij (1^) 

and a second term similar to the dynamic Smagorinsky models such as ([5]). How- 
ever, in the case of ([T^ . the term proportional to Lij derives from the Germano 
identity directly. It is well-known from the literature, that similarity models better 
reproduce the SG tensor components, whereas models based on the concept of eddy 
viscosity ((T|), such as dynamic Smagorinsky model, better reproduce dissipated en- 
ergy [5]. Tensor may be able to better reproduce both SG tensor components 
and dissipated energy. 
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